For  Reference 


not  to  be  taken  from  this  room 


The  University  of  Alberta 
Printing  Department 
Edmonton,  Alberta 


UNIVERSITY  OF  ALBERTA 
LIBRARY 


Regulations  Regarding  Theses  and  Dissertations 


Typescript  copies  of  theses  and  dissertations  for  Master’s  and  Doctor's 
degrees  deposited  in  the  University  of  Alberta  Library,  as  the  official  Copy  of 
the  Faculty  of  Graduate  Studies,  may  be  consulted  in  the  Reference  Reading  Room 
only. . 


A  second  copy  is  on  deposit  in  the  Department  under  whose  supervision  the 
work  was  done.  Some  Departments  are  willing  to  loan  their  copy  to  libraries, 
through  the  inter-library  loan  service  of  the  University  of  Alberta  Library. 

These  theses  and  dissertations  are  to  be  used  only  with  due  regard  to  the 
rights  of  the  author.  Written  permission  of  the  author  and  of  the  Department 
must  be  obtained  through  the  University  of  Alberta  Library  when  extended  passages 
are  copied.  When  permission  has  been  granted,  acknowledgement  must  appear  in  the 
published  work. 

This  thesis  or  dissertation  has  been  used  in  accordance  with  the  above 
regulations  by  the  persons  listed  below.  The  borrowing  library  is  obligated  to 
secure  the  signature  of  each  user. 


THE  UNIVERSITY  OF  ALBERTA 


EFFECTS  OF  THE  INDUCED  MAGNETIC  FIELD 
ON  THE  PLASMA  COLUMN 


by 

DAVID  SEALE 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES 
IN  PARTIAL  PULFILMENT  OF  THE  REQUIREMENTS 
FOR  THE  DEGREE  OF 
MASTER  OF  SCIENCE 

DEPARTMENT  OF  PHYSICS 


EDMONTON,  ALBERTA 


AUGUST,  1967 


Digitized  by  the  Internet  Archive 
in  2019  with  funding  from 
University  of  Alberta  Libraries 


https  ://arch  ive  .org/detai  Is/Seale  1 967 


UNIVERSITY  OF  ALBERTA 
FACULTY  OF  GRADUATE  STUDIES 

The  undersigned  certify  that  they  have  read,  and 
recommend  to  the  Faculty,  of  Graduate  Studies  for  acceptance 
a  thesis  entitled  Effects  of  the  Induced  Magnetic  Field  on 
the  Plasma  Column,  submitted  by  David  Seale,  in  partial  ful 
filment  of  the  requirements  for  the  degree  of  Master  of 


Science . 


— 


ABSTRACT 


Beginning  with  Euler's  equation  of  hydrodynamics 
and  the  equation  of  continuity,  using  the  method  set  forth 
by  Johnson  in  1962,  the  steady  state  profile  and  the  solu¬ 
tion  of  the  helical  mode  are  found  for  the  plasma  column 
in  parallel  electric  and  magnetic  fields.  In  this  calcul¬ 
ation,  linear  terms  in  the  induced  azimuthal  magnetic  field 
are  retained.  In  previous  work,  this  induced  field  was 
neglected  altogether. 

In  the  steady  state,  the  well  known  pinching  due 
to  the  azimuthal  field  is  observed.  In  the  helical  solution 
the  parameters  of  the  helix  are  calculated  as  a  function  of 
the  critical  longitudinal  magnetic  field  at  which  the  helical 
state  is  stable  and  it  is  found  that  inclusion  of  the  azi¬ 
muthal  field  terms  results  in  a  deviation  from  Johnson's 
theory  for  small  values  of  the  critical  longitudinal  mag¬ 
netic  field,  when  the  latter  does  not  completely  overshadow 
the  azimuthal  field. 

Experimental  observation  of  this  deviation  would 
require  experiments  especially  designed  for  this  purpose--. 
Under  normal  conditions  of  current  and  density,  the  effect 
of  the  azimuthal  field  is  negligeable .  To  date  no  such 
experiment  has  been  performed. 
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NOTATION 


ion-electron  density 

flux 

temperature 

frictional  relaxation  time 

ion-electron  creation  rate 

particle  mass 
magnetic  field 

electric  field 

=  v  ,  -  v  =  the  sum  of  the  absolute 

z+  z  - 

values  of  the  longitudinal  velocities.  vz 
is  positive  for  positive  electric  field  and 
is  proportional  to  the  current. 

average  ion-electron  density 

absolute  value  of  the  electronic  charge 

velocity  of  light 


mobility 


7o'  71J 

successive  zeros  of  J 

o 

V  v  ••• 

successive  zeros  of 

R 

radius  of  glass  tube 

D 

coefficient  of  diffusion 

parameters  associated  with  the  azimuthal 

magnetic  field 

0) 

frequency  of  the  helix 

K 

wave  number  of  the  helix 

+  and  -  refer  to  ions  and  electrons.  The  cgs-esu  system 
of  units  is  used  throughout. 
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Chapter  1;  Introduction 


1.1  Description  of  the  system 

The  system  to  be  studied  consists  of  a  plasma., 
that  is  a  partially  ionised  gas,  in  parallel  electric  and 
magnetic  fields.  The  plasma  is  contained  in  a  cylindrical 
glass  tube  whose  diameter  is  of  the  order  of  a  few  centi¬ 
meters  and  whose  length  is  of  the  order  of  a  meter.  The 
effect  of  the  ends  is  neglected.  The  electric  field  is 
produced  by  electrodes  at  the  ends  of  the  tube  and  the  mag¬ 
netic  field  by  a  coil  wound  around  the  tube.  The  plasma 
state  is  obtained  in  a  sustained  electric  discharge  in  the 
tube.  This  discharge  is  produced  by  the  longitudinal  el¬ 
ectric  field  which  also  serves  to  vary  the  current  in  the 
discharge.  The  gas  pressure  is  of  the  order  of  1  mm  Hg 
and  the  gas  density  is  approximately  1C>  molecules/cmr  . 

In  the  range  of  electric  fields  used„  the  charged  particle 

IS  ,  s 

density  is  of  the  order  of  10  parti cles/crrm . 

The  mobilities  of  the  ions  and  electrons  is  con¬ 
sidered  to  be  determined  by  collisions  with  neutrals  only 
since  neutrals  outnumber  charged  particles  by  several  orders 
of  magnitude.  Since  the  neutral  gas  pressure  is  uniform 
throughout  the  tube.,  the  ionic  and  electronic  mobilities  are 
regarded  as  constant  with  respect  to  all  variables.  These 
mobilities  characterize  the  gas  under  consideration  in  the 
calculation.  As  an  example.,  the  mobilities  of  ionized  helium 


' 

. 
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are  used  for  the  numerical  work. 

The  walls  of  the  tube  act  as  an  ion-electron  pair 
sink.  This  is  reflected  in  the  boundary  condition  that  the 
charged  particle  density  at  the  walls  is  zero.  This  results 
in  a  negative  radial  gradient  of  charge  density.  This  leads 
to  an  outward  radial  diffusion  of  ions  and  electrons  whose 
interaction  with  the  longitudinal  magnetic  field  produces  an 
azimuthal  force  on  the  particles.  The  resulting  azimuthal 
motion  and  the  longitudinal  motion  of  the  charged  particles, 
due  to  the  external  electric  field,,  combine  to  produce  a 
helical  motion.  This  helical  state  is  studied  in  chapter  3. 

To  avoid  a  buildup  of  charges ,  the  radial  fluxes  of 
ions  and  electrons  must  be  equal.  Due  to  the  difference  in 
the  mobilities  of  ions  and  electrons.,  we  must  have  an  internal 
radial  electric  field  to  draw  ions  to  the  walls  and  retard  the 
electrons.  The  densities  of  ions  and  electrons  are  assumed 
equal  at  every  point.  The  possibility  of  an  internal  electric 
field,  however,  is  not  excluded.  These  are  the  assumptions 
of  quasi-neutrality,  which  is  justified  in  chapter  2. 

The  neutral  gas  remains  at  approximately  room  tem¬ 
perature.  However  the  temperature  of  the  electrons  is  sev¬ 
eral  orders  of  magnitude  greater  in  the  range  of  pressures 
and  electric  field  used.  Due  to  the  infrequency  of  collisions 
between  electrons  and  neutrals,  the  electrons  acquire  much 
energy  from  the  electric  field  before  transferring  energy  to 
the  gas  in  a  collision.  The  average  energy  so  gained  must 
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equal  the  average  loss  at  each  collision.  Since  most  col¬ 
lisions  are  elastic.,  the  transferred  energy  is  a  small  frac¬ 
tion  of  the  total  energy  (Howatson  1965) •  Therefore  the  av¬ 
erage  energy  of  the  electrons  is  great  in  an  electric  fie  2d 
and  consequently  their  temperature  is  high.  The  ions,  on  the 
other  hand,  experience  frequent  collisions  and  exchange  ap¬ 
preciable  kinetic  energy  with  neutrals.  Therefore  their  temp¬ 
erature  remains  close  to  the  temperature  of  the  gas.  Since 
the  charged  particle  temperatures  are  determined  by  the  freq¬ 
uency  of  collisions  with  neutrals  and  since  the  density  of 
neutrals  is  constant  in  space,  we  assume  the  charged  particle 
temperatures  independent  of  spatial  coordinates. 

Due  to  the  sink  at  the  walls,  a  concentration  of 
charges  is  observed  in  the  axial  region.  This  effect  corres¬ 
ponds  to  the  steady  state  solution  given  in  chapter  2.  As 
the  externally  applied  longitudinal  magnetic  field  is  increas¬ 
ed,  a  screw-type  state  is  observed.  This  helical  state  is 
observed  to  be  unstable  for  sufficiently  large  magnetic  fields. 
The  mathematical  description  of  this  helical  state  is  given 
in  chapter  3. 

The  variables  of  this  problem  are  the  density  of 
charged  particles,  the  magnetic  field  and  the  electric  field. 
The  longitudinal  electric  and  magnetic  fields,  since  they  are 
externally  applied,  are  treated  as  independent  variables. 

The  internal  radial  electric  field,  the  density,  and  the  in¬ 
duced  magnetic  fields  are  the  dependent  variables. 
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1.2  Previous  work  in  this  field 

The  existence  of  the  critical  magnetic  field  was 
first  described  by  Guthrie  and  Wakerling  in  1949*  The  first 
satisfactory  theoretical  treatment  of  the  problem  was  per¬ 
formed  by  Kadomtsev  and  Nedospasov  in  i960.  Their  calcul¬ 
ation  was  refined  by  Johnson  and  Jerde  in  1962  and  their 
treatment  remains  the  most  complete  on  the  subject  to  date. 

The  theory  of  Kadomtsev  and  Nedospasov  begins  with 
the  equation  of  continuity  and  the  equations  of  motion,  and 
by  dictating  a  perturbed  radial  dependence  for  the  electron 
density  and  potential,  using  the  zeroth  order  expressions, 
obtains  a  dispersion  relation.  This  is  done  by  a  normal 
mode  analysis. 

Johnson  generalizes  this  method  by  removing  the 
restriction  of  having  to  specify  the  radial  dependences  of 
density  and  potential,  although  in  his  actual  calculations, 
he  does  not  go  beyond  the  zeroth  order.  The  advantage  of 
his  treatment  is  that  it  could  be  easily  be  adapted  to  in¬ 
clude  higher  orders. 

Furthermore,  Kadomtsev  and  Nedospasov  neglect  ionic 
diffusion  and  interaction  with  the  magnetic  field,  thereby 
limiting  the  application  of  the  theory  to  cases  of  large 
electron-ion  mobility  ratios.  Johnson,  however,  makes  no 
such  approximations,  and  his  calculation  remains  symmetric 
with  respect  to  ions  and  electrons. 

None  of  the  theoretical  work  in  this  field  has  ever 
included  the  term  due  to  the  internal  magnetic  field  in  the 
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azimuthal  direction.  The  explanation  for  this  is  that  first  * 
the  longitudinal  magnetic  field  normally  overshadows  the  azi¬ 
muthal  field  and  second,  inclusion  of  this  term  complicates 
the  algebra  enormously. 

As  a  consequence  of  the  omission  of  this  term  we  can 
expect  that  the  results  of  the  above  authors  would  be  unreal¬ 
istic  for  small  critical  magnetic  fields  when  the  azimuthal 
field  becomes  comparable  to  the  longitudinal  field. 

1.3  The  improvement  to  be  made 

The  object  of  the  present  calculation  is  to  include 
the  term  due  to  the  induced  azimuthal  magnetic  field.  John¬ 
son's  method  is  extended  to  include  this  term  and,  as  it  will 
be  shown  later,  the  results  for  small  critical  magnetic 
fields  are  rendered  physically  realistic. 

The  additional  term  is  assumed  to  be  small,  and 
only  the  linear  terms  are  kept.  The  generality  of  Johnson's 
work  is  forsaken  in  taking  only  zeroth  order  terms  from  the 
start  because  of  the  complexity  of  the  algebra.  However, 
since  Johnson  has  carried  out  numerical  calculations  only 
to  the  zeroth  order  anyway,  the  results  of  this  calculation 
reduce  to  his  if  the  parameter  associated  with  the  azimuth¬ 
al  magnetic  field  is  reduced  to  zero.  The  difference  is 
that  the  present  calculation  could  not  easily  be  adapted 
to  include  higher  orders. 

In  the  next  chapter,  orders  of  magnitude  are  given 


, 
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for  the  quantities  involved  and  it  is  found  that  the  effect 
of  the  additional  term  is  not  negligeable  under  certain 
circumstances.  The  additional  term  is  most  important  in 
cases  of  large  current  and  of  large  average  density.  Since 
Johnson's  theory  requires  no  upper  bound  on  these  quanti¬ 
ties,  there  is  a  region  where  the  results  of  this  calcul¬ 
ation  may  be  observed.  The  system  may  be  placed  in  cond¬ 
itions  where  the  effect  of  the  azimuthal  field  may  be  ob¬ 
served  without  sacr  ificing  the  validity  of  Johnson's  basic 
theory . 

1.4  Outline  of  the  method 

The  treatment  is  divided  into  two  parts.  First 
the  steady  state  is  studied.  The  steady  state  is  defined 
by  the  condition  that  the  density  at  every  point  in  space 
remains  constant  in  time.  Beginning  with  the  equations  of 
motion,,  which  are  derived  from  Euler's  hydrodynamic  equa¬ 
tion,  and  the  continuity  equation,  a  solution  in  the  form 
of  a  series  of  Bessel  functions  of  order  zero  is  obtained. 
The  coefficients  of  this  series,  except  for  the  first,  are 
zero  for  zero  azimuthal  magnetic  field,  yielding  Johnson's 
solution,  which  is  a  single  zero-order  Bessel  function. 

The  coefficient  of  the  first  order  term  is  calculated. 

In  the  second  part,  beginning  with  the  same  equa¬ 
tions,  the  helical  mode  is  studied.  A  solution  in  the  form 
of  a  series  of  Bessel  functions  of  order  one  is  substituted 


. 


7. 


and  from  the  relation  between  the  coefficients ,  using  a 
stability  criterion,  an  equation  containing  the  wave¬ 
length  of  the  helix,  the  critical  longitudinal  magnetic 
field,  and  the  longitudinal  electric  field  is  obtained. 
The  equation  also  contains  the  azimuthal  magnetic  field 
as  a  parameter  and  reduces  to  Johnson's  corresponding  eq¬ 
uation  if  the  azimuthal  field  is  reduced  to  zero. 
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Chapter  2:  Basic  Equations  and  steady  state  solution 
2.1  Introduction  of  the  basic  equations 

We  obtain  the  equations  of  motion  of  the  system 
from  Euler's  equation  of  hydrodynamics.  Euler's  equation 
may  be  written 

p  {  +  (v-grad) ]  v  +  V  p  =  P 

where  p  is  the  mass  density,  v5”  the  velocity,  p  the  pressure 
— > 

and  F  the  external  forces  (Sommerfeld) .  In  the  further 
use  of  this  equation,  we  neglect  (v*grad)v  in  comparison 

with  grad  p  or  the  external  forces. 

\ 

We  treat  the  plasma  as  two  fluids  and  write  an 
Euler  equation  for  each  of  the  fluids,  electronic  and  ionic. 

We  may  write  these  equations 

— > 

v±  * 
n±m±  —  +  P±  =  F± 

where  n+  is  the  particle  density  and  m+  the  particle  mass. 
The  4-  and  -  indice-s  refer  to  ions  and  electrons. 

The  external  forces  acting  on  these  fluids  are  the 
frictional  force  and  the  forces  due  to  the  electric  and 
magnetic  fields.  These  may  be  written 

n  +m  +v  + 

^friction±  ~~  t+ 
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— >  -4 


p  _ 

electric± 

±  n+  e  E 

— > 

.  — ^ 

F  = 

magnetic  ± 

±  n  ,  e  v 
+ 

where  t+  is  the  relaxation  time  responsible  for  the  fric- 

— >  — ^ 

tional  force*  E  the  electric  field*  B  the  magnetic  field* 
and  e  the  absolute  value  of  the  electronic  charge.  Though 
the  magnetic  force  is  quadratic  in  the  longitudinal  velocity 
component  v  *  this  term  is  larger  than  the  term  involving 

Z 

v* (grad  v)  since  the  latter  must  involve  either  the  radial 
or  azimuthal  velocity  which  is  much  smaller  than  v  . 

We  note  that  the  azimuthal  magnetic  field  interacts 
with  the  longitudinal  current  to  create  a  radial  force . 
Therefore  the  azimuthal  field  will  modify  the  radial  motion 
whose  interaction  with  the  longitudinal  magnetic  field  re¬ 
sults  in  the  azimuthal  motion*  as  seen  in  chapter  1.  Con¬ 
sequently*  the  azimuthal  field  will  affect  the  helix. 

To  solve  the  problem*  we  make  the  assumption  of 
quasi -neutrality .  The  densities  of  ions  and  electrons  are 
regarded  as  equal  everywhere  without  excluding  the  possib¬ 
ility  of  an  internal  electric  field.  This  assumption  is 
based  on  the  fact  that  departures  from  neutrality*  insig¬ 
nificant  from  the  point  of  view  of  density*  create  strong 
electric  fields.  This  assumption  permits  the  removal  of 
the  +  and  -  indices  from  n. 
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We  assume  both  ionic  and  electronic  gases  ideal 
and  write  the  equation  of  state  as 


p+  =  n  kT+ 


A  boundary  condition  of  the  problem  is  dictated  by 
the  sink  at  the  walls.  This  is  written  as 


n(R) 


0 


where  R  is  the  radius  of  the  tube. 

The  other  hydrodynamic  equation  which  we  must  use 
is  the  continuity  equation  for  a  gas  whose  particles  are 
not  conserved,  which  may  be  written  as 


dn  .  r-,  t  — >\ 

St  +  v'(n  v) 


=  n| 


where  is  the  number  of  pairs  created  per  electron  per  unit 
time.  Due  to  the  dominance  of  recombinations  at  the  boundary, 
we  neglect  recombination  and  assume  only  regeneration  in  the 
volume  of  the  tube.  The  rate  of  production  of  ion-electron 
pairs  is  due  to  ionization  of  the  neutral  molecules  by  col¬ 
lision  with  electrons  and  is  directly  proprtional  to  the  num- 
ber  of  electrons.  Therefore,  the  rate  of  regeneration/ cm 
is  n£,  where  n  is  the  number  of  electrons. 

Since  the  frequency  of  the  helix  is  much  less  than 
1/t,  we  neglect  the  time  derivative  of  the  velocity  in  Euler's 
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equation  (Spitzer,  Johnson) .  This  assumption  is  equivalent 
to  assuming  the  system  in  equilibrium  with  the  external 
forces  at  every  instant. 

The  basic  equations  may  then  be  written 


f-L.  —  P+  ( 


kT 


Vn±nE±  T,  X  B  ) 


a) 


!£  +  v-r± 


=  n  | 


(2) 


p+  =  n  kT  + 


(3) 


where  we  define  T+  and  as  follows 


— >  — » 

T+  =  nv+ 


ex 


+ 


+ 


m 


+ 


We  note  that  the  continuity  equation  requires  that 


V-F+  =  v-r 


' 
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2.2  Steady  state  solution 

We  define  the  steady  state  by  the  following  con¬ 


dition 


dn 

3t  = 


0 


Writing  the  r  components  of  the  equations  of  motion,  we  ob¬ 
tain 


1 

-+1-. 

-n  E 

-  n 

(v  B  _  -  v  B  ) 

'  -cp  oz  -z  oep' 

kT_ 

dn  | 

~5r  _  =  n  v-r 

or 

e 

-Hi+l 

n  E 

+  n 

(v.  B  -  v,  B  ) 

v  +<p  oz  +z  ocp7 

kT 

1?  =  n  v+r 

(4) 

or 

e 

which  are  first  order,  partial  differential  equations  in  n. 
One  of  the  boundary  conditions  for  equations  (4)  is  dictat¬ 
ed  by  the  sink  at  the  surface.  That  is 

n(R)  =  0 

Because  of  the  cylindrical  symmetry  in  the  steady 
state,  to  avoid  charge  separation,  we  must  have 


Solving  for  v,  and  v 

°  +cp  -ep 


equation  (2)  we  get 


from  the  cp  component  of 


v 

+<P 


\L  ,  V 

+  r 


^-vr 


(5) 


' 
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since  E  and  -^3- 
ocp  dcp 

B  must  be  zero 
or 


must  be  zero  for  cylindrical  symmetry  and 
since  it  must  satisfy  Maxwell's  equation 


V  •  B  =  0 


and  also  be 
We 

sion  of  B 

ocp 


where  n  is 


cylindrically  symmetric. 

use  the  lowest  term  in  the  power  series  expan- 
,  that  is  the  term  linear  in  r.  We  assume 
n  e  v  r 


B  ~  — 

ocp 

the  average 


- § -  (6) 

C 

ion-electron  density  and 


v 


z 


V 


z  - 


which  is  proportional  to  the  current. 

Between  equations  (4)  we  eliminate  EQr.  This  yields,, 
using  substitutions  (5)  and  (6),  and  solving  for 


r 


r 


2  - 

v  n  e  r 
z 


2 

c 


2  kT  ^n  -I 
e  "5r  J 


(7) 


where 


M-  = 


1 


( M-  ■  +  M-  )  t  1 


■+  '  +  B  ) 

\i+\i  ozJ 


and 


T  +  T_ 

T  =  -  =  average  temperature 


. 

■0.9 
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Substituting  the  expression  for  from  equation 
(7)  into  the  continuity  equation  (2)  we  obtain  the  second 
order  differential  equation 


d2n 

8r^ 


+ 


[  1  + 
r 


2  -  l 
v  n  e 
z _ 

2c2  kT 


)] 


8n 

■5? 


+  [ 


2  - 

v  n 
z 


+ 


kT 


e  I 

2p"kT 


] 


n 


0 


Since  n  depends  only  on  r  we  may  replace  the  partial  deriv¬ 
atives  by  total  derivatives. 

We  define  the  dimensionless  variables  p  and  N  as 

follows 

p  =  r/R  N  =  n/n 

where  R  is  the  radius  of  the  tube  and  n  is  the  average  el¬ 
ectron-ion  density.  The  differential  equation  may  then  be 
written 


d2N 

dj7 


+ 


(^+i£-)^+(A+B)N=0 


dp 


(8) 


where 


A  = 


2  -  2 
v  n  e  R 
z 

2 

c  kT 


B  = 


e  R2? 

2jjL"  kT 


and  N(R)  =  0 

We  note  that  A  and  B  are  dimensionless. 

If  we  put  A  equal  to  zero,  which  corresponds  to 
putting  B  to  zero,  equation  (8)  reduces  to  an  exact  Bes¬ 
sel  equation,  whose  solution  is  a  Bessel  function  of  order 
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zero  and  whose  first  zero  occurs  at  R,  which  corresponds 
to  Johnson's  solution  for  the  steady  state.  We  can  ex¬ 
amine  the  effect  of  BQCp  by  writing  the  solution  as  a  ser¬ 
ies  of  zero  order  Bessel  functions  whose  first,  second, 
etc.,  zeros  occur  at  r  =  R. 

The  solution  is  written  as 


N  =  N  [  J  (7  p)  +  an J  (y_p)  +  ...  ] 
o  o  wor  lowlr/ 


where  7,7,,  ...  are  the  successive  zeros  of  J  . 
o  1  o 

Substituting  this  solution  into  equation  (8), 
multiplying  by  J^y^p)  p  dp  and  integrating  over  p  from 
zero  to  one,  we  obtain  an  expression  for  a^  whose  linear 
term  in  A  is 

1  2 

A  /  p  joKp>  dp 

1 

2  (  ~  )  /  Jq(7iP)  P  dP 

The  numerical  values  of  the  integrals  give 

a2  =  0.0346  A  +  . . . 


Therefore  the  lowest  order  expression  for  a^  becomes 


-  2 

n  e 


~  0.0346 


2 

c  kT 


and  the  solution  may  be  written 


v2  -  02  r2 

N  -  N0  [  J0(70p)  +  0.0346  -  -VUP)  ]  (9) 

c2  kT 


The  orders  of  magnitude  of  the  quantities  in  the 
above  solution  are 


v  ~ 

z 

n  ~ 

e  ~ 

R  ~ 

C  'v 

k  ~ 

T 


,6 


10  cm/sec 
10^  /cm  ^ 


in'10 

10  esu 


1  cm 


1010  cm/sec 
- 1 6 

10  ergs/degree  K 


4 

10  degrees  K 


(Brown  1959) 


which  give  an  order  of  magnitude  of  10  1  for  A.  For  the 
approximation  of  taking  only  the  linear  term  in  the  ex¬ 
pansion  for  N4  we  must  have  a^  «  1.  That  is  we  must  have 
A  «  30.  To  plot  the  profile  so  that  the  effect  of  the 
azimuthal  field  is  not  negligible,  we  take  A  =  1.  In 
this  case  the  solution  is 


N  =  Nq  [  J0(70P)  +  0.0346  J0(7XP) ] 


(10) 
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where  the  numerical  values  of  the  7's  are 
7o  =  2.405 
7X  =  5.520 

For  this  example  the  results  of  the  calculation  are  shown 

in  figure  1.  From  the  graphs  we  observe  that  the  calculated 

result  is  reasonable  since  it  is  known  that  the  effect  of 

B  in  semiconductors  is  to  pinch  the  column, 
ocp 


P 

figure  1 

Graph  of  equation  10  (solid  curve)  for 
parameter  A  =  1,  and  Johnson's  solution, 
l(70p)  (dashed  curve) 
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Chapter  3:  Helical  solution 

It  has  been  established  both  theoretically  and 
experimentally  that  a  helical  mode  exists  for  this  system. 
Johnson  has  developed  the  theory  and  performed  calcula¬ 
tions  for  the  helical  mode.  In  this  section,  we  shall 
examine  the  effect  produced  by  the  terms  due  to  Bo(^  which 
Johnson  has  neglected. 

For  longitudinal  magnetic  fields  below  a  critical 
value,  which  we  shall  calculate.,  the  helical  mode  decays. 
Above  the  critical  field  the  helical  state  is  unstable. 

The  imaginary  part  of  the  frequency  depends  critically 
on  the  longitudinal  magnetic  field. 

We  begin  with  the  same  basic  equations  as  in  the 
previous  chapter.  The  equations  of  motion  are  as  before 

r ±  =  q±[  ±  n  E  ±  ?±  X  B  -\7(kT±  n)  ]  (10) 

and  the  equation  of  continuity  is 

V-?±  =  -  |S-  +  ni  (11) 

First  we  must  solve  equation  (10)  for  an  explicit 
— ) 

expression  for  T+  which  may  be  substituted  into  equation 
(ll) .  To  do  this  we  write  equation  (10)  as  three  equations 
in  the  cylindrical  components  of  T+.  From  the  expressions 
of  its  components  we  can  construct  a  vectorial  expression 

— ^ 

for  r+.  This  expression  may  be  written 


. 

,r . 

. 
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— > 
r 


|i.9  n  B*E  B 


kT 

O  ^  ^  ^  ^  ^ 

q  +  -— 1  B  B-'Vn  ±  n  [  E  +  B  X  E  q 


+ 


+ 


2  2 
1  +  m^+b 


ki+kT  + 


[  q,B  X  Vn  +  Vn  ] 


4* 


2  2 

1  +  q+B 


(12) 


Combining  equations  (ll)  and  (12)  to  eliminate  T+  we  obtain 


bn 

3t 


+  n£  =  V’{[±q+  nB-E  B  - 


M«1  kT  + 


->  — > 


->  - 


B  B-Vn  ±  n(E  + 


_>  q+kT  + 

q+B  X  E)  ±  (q+B  X  Vn)  +  Vn  ]  / 


(1  +  q+B2)  } 


(13) 


This  equation  may  be  simplified  by  using  Maxwell's  equation 


V  •  B  =  0 


and  the  identity  from  vectorial  analysis 

i 

V  X  V  n  =  0 


Since  the  azimuthal  field  is  normally  much  smaller  than 

2  2 

the  external  longitudinal  field,  we  replace  1  +  q+B  by 


20. 


1  +  11+  B  .  Using  these  conditions,  we  obtain  from  equa- 

—  U  kj 

tion  (13) 


(1  +  4boZ)(-  ||  +  n0 


BaVn  B*E  +  B*V(B*E) 


\±l  kT±  B-V(B°Vn)  ±  |i+E*Vn  -  [±± 


n 


E?-  (V  X  B)  -  nEh  (V  X  E?) 


— > 


+  ( B  X  E )  •  Vn 


2 


2 


±  q1!  kT+  Vn  •  ( V  X  B)  -  p±  kT+  V  n 


To  treat  this  problem,  we  define  a  new  variable  np  as  fol¬ 
lows 

n-i  =  n  -  n 
1  o 

where  n  is  defined  as 
o 

n  =  N  n 
o 

where  n  is  the  average  density  and  N  is  the  steady  state 
solution  given  in  section  2.2.  We  assume  that  n^  is  small 
and  depends  on  the  coordinates  and  time.  Also  we  must  put 


— >  — >  -» 

E  =  E  +  E, 
o  1 

— >  — > 

where  Eq  is  the  electric  field  in  the  steady  state  and  E^ 

is  time  and  coordinate  dependent  and  assumed  small.  Re¬ 
placing  n  and  E  with  the  above  expressions  in  (14),  from 
which  has  been  subtracted  the  steady  state  equation,  we  ob¬ 
tain 


p  p  dn-.  P  P  clE-, 

(1  +  U±B^z)(  -  ^  +  )  =  ±H?  n„Bf 


dn. 


‘13  v 

E  v;  ±  |i3n  BE  —  -5— 
oz  d^z  o  oz  oz  r  dr 


^  4  (r  B0(p)  - 


±  'o  oz  oz 
kT 
±  ~e 


3  2 

±  ii-tB  * 
oz 


±  2  (l  d  , 

®oz  \r  dr^rdz  ^ 


■ 
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-n2 

a  n 

J_  JT.  -  -  J_ 

r  dcpdz 


>2 
d  n 


1 


a?  /  "±n° 


r  Br  (rEi + 


Jlr'  r  Bcp 


as,  dE, 

lep  +  lz 


oz 


^±no 


B  d  ( rE ,  „ )  ,  aE, 

oz  v  lcp y  1_  lr 

r  Br  ""  oz  r  dcp 


+ 


+ 


an, 

1  T? 

dz 

oz 

r  dr 

pi±kT± 

i  s2ni 

OCp 


E 


oz 


2 

pi+  n 


o 


r  Br  (rn]_Eor) 
r  dcp  Boz  Eor^±  ^ 


-M- j.  —  B 


'  ±  r 


Brd'cp 


2_  a2(rn1) 

oz  '  B±r2  BrBep  Boz 


V2n- 


pi+kT  + 


(15) 


We  now  assume  a  helical  form  for  the  density  n 


1 


and  the  potential  V-^  of  the  field  E^ .  This  is  expressed  as 


n^  =  f(r)  exp  [  i(Kz  +  mcp  +  cot)] 


V  =  g(r)  exp  [  i(Kz  +  mcp  +  cot)] 


Since  we  wish  to  study  the  simple  helix.,  we  set  m  =  1. 

The  electric  field  may  be  expressed  in  terms  of  the  poten¬ 
tial  according  to  the  following  set  of  equations: 


E,  =  - 
lr 

ag  ^i 

Br  g 

Elq»  “  *  1 

vi 

r 

E  = 
lz 

-  iKV1 

aE, 

lr  = 

a2g 

ar2 

a(rEld  = 

-  i 

dg 

dr 

dE, 
lz  _ 

dr 

- 

dr 

dr 

dr 

1  lr 

i  ag 

,  dE, 

1  lep 

g 

9Elz 

=  Kg 

r  dcp 

r  dr 

r  Bcp 

2 

r 

dcp 

aE, 

lr 

-  ¥■ iK 
dr 

dE, 

lep 

Kg 

dE, 

lz 

=  K2g 

dz 

dz 

r 

dz 

' 
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Defining  the  quantities  |i,J.  and  D'+  as  follows 


M-i  = 


D 


M*± 

- 

1  +  q+B 

oz 


D± 

2~2~ 

1  +  Moz 


Li+kT  + 


- S-? 

1  +  ^±Boz 


where  D+  is  the  coefficient  of  diffusion  (see  Einstein's 
relation  between  the  mobility  and  the  diffusion  coeffic¬ 
ient  ) , replacing  B  by  its  lowest  order  expression,,  that 


is  the  term  linear  in  r 

B 


n  vz  e  r 


ocp 


and  defining 


B 

.  ocp  n  v  e 

r  2 


from  (15)  we  obtain 

t  1  +  n  *  I  4-  (rE  f)  -  [(io)  -  i)  ±  |i+F  IK  +  D,K2 

±  r  dr^  dr7  r>  r\-r>  >  LV  s  ^  ±  OZ  ± 


D 


r  dr  '  or 


+  D'/r 


B 


oz  . 


iE  LLuii  ]  f  +  [rF 


or'  '±^± 


oz  dr 


df  HUH  +  M-i  41BozEozf 


2, 


-  2KB^r,[i+D|f  +  2(E  U±  +  D±iK) i-i  +f  -  iKrE  Ji^±f  ]  $ 


oz' 


n 


o 


+  ^±?  3F(rnoil)  +  [±^±  y 


+  [  Hi  4  2n  B  K  +|a;n±noiK2+  ni|X±i  )Kr  ]  g  $  (17) 


or 

Q  dn  B 

1  H±noKd  -  IX in ±i(^)  -4-  ]  g 


dn 


Since  all  variables  depend  only  on  v,  we  have  replaced 
the  partial  derivatives  by  total  derivatives. 
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In  equation  ( 1 7 )  terms  linear  in  0  have  been  kept 
whereas  term^ containing  higher  powers  of  $  have  been  ne¬ 
glected.  We  now  define  a  new  variable  I  such  that 


g(r) 


Air) 

Vr) 


(18) 


Since  g(R),  which  is  essentially  the  amplitude  of  the  pot¬ 
ential  of  the  electric  field  at  the  boundary,  must  be  fin¬ 
ite  and  since  n  (R)  =  0  due  to  the  ion-electron  sink  at  the 
boundary,  I(R)  must  be  zero.  Rewriting  equation  (17)  in 
terms  of  I  and  f  we  obtain 


Di  i  4r  (r  g  )  +  n;  i  fj;  (rEorf)  -  [(ia>  -  |)  ±  IK 


±  r  dr 


'±  r  dr 


p  p  B 

+  D±K^  +  D+/ r  -  -|5  iEor  p.±p.j.  ]  f  +  [  rE 


df 
oz  dr 


±  oz 


M.+M4  + 


f  -  2KB  p,  Di  f  +  2  (E  p,  +  D.iK)  p,f 

OZ  ±  ±  v  OZ.  ±  ±  '  ± 


iB  E 

± ±  oz  oz 

i  f  ]  O  =  +  \±\  ± 

+  K  ±  r  dr  v  dr 


KriE  ijl  .  |jl  1 
or  ±  ± 


1  d  / 

(n 


dl  )  ±  ,  1  d(lrp) 

'  r  dr 


±  7  ±  n+  K2  -  np±  225.  ]  n  +  [  ±  uin2  2BqzK  +  2d>±lK 

r 


+  [ 

+  ipKr  ]  <&& 


where  p(r) 


1_ 

n 

o 


dn 
_ o 

dr 


We  know  from  Johnson's  work  that  the  approximate 
solution  to  this  problem  without  the  terms  in  0,  is  a  Bes- 


. 


' 
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sel  function  of  order  one,  whose  first  zero  occurs  at  the 
boundary.  Since  we  can  vary  the  parameter  0  over  a  wide 
range  without  losing  the  validity  of  Johnson's  theory, 
we  note  that  the  introduction  of  $  is  consistent  with 
Johnson's  approximations.  We  may  write  the  solution  to 
the  problem  dealing  with  the  terms  in  0  as  the  sum  of 
first  order  Bessel  functions  whose  successive  zeros  oc¬ 
cur  at  r  =  R.  We  write 

f  =  f0  Ji(Kor)  +  fl  Ji(Klr)  +  • • • 

(21) 

l  =  lQ  q(Kr)  +  •«!  +  .  .  . 

where  K  ,  K,  ,  ...  are  the  first,  second,  . . .  zeros  of  J,  . 
o  1  1 

Substituting  expressions  (21)  into  equation  (19)*  multi¬ 
plying  by 

J_  (K  r)  r  dr 
1 v  o  ' 

D± 

and  integrating  over  r  from  r  =  0  to  r  =  R,  we  obtain  two 

equations  relating  fQ  and  l ^  with  terms  in  f^  ...  and 

...  .  All  terms  containing  f^,  ...  and  ...  are 

multiplied  by  integrals  of  the  form 

R  R 

/  J1(KQr)  q(qr)  r  dr  ,  /  q(KQr)  q(qr)  r  dr,  .. 

If  the  range  of  integration  were  infinite,  all  these  int¬ 
egrals  would  vanish,  a,nd  this  method  would  reduce  to  the 
standard  transform  method.  However  in  the  case  of  a  fin- 


ite  range  of  integration,  these  integrals  are  finite  but 
much  smaller  than  the  integral 
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0 


/  J^(KQr)  r  dr 


which  occurs  in  the  coefficients  of  the  f  and  &  terms. 

o  o 

This  is  due  to  the  oscillation  of  the  integrands  which 

contain  Bessel  functions  of  different  arguments.  Also, 

since  f  *  . . .  and  i^,  . . .  are  only  corrections,  they 

must  be  small  compared  to  f  and  i  .  Consequently,  we 

neglect  all  terms  containing  f^,  ...  and  ,  ...  and  are 

left  with  two  equations  in  f  and  i  of  the  form 

00 


A  f  +  B  i 
o 


o 


0 


where 


o  ±  3 


+ 


iK  +  — rK  kT 


2KB  \i2,  kT, 
oz  ±  ± 


2Eozn±  ±  2i|x±  KkT±  +  i  p.±  K  IgCt  ] 


and 


B 


±  [  K2  +  IQ  + 
o  3 


+  2|jl±  iK  +  q±iK  Ip  ] 


■ 
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and  where 


a  = 


M'x  kT  -  |a '  kT 

+  U~ _ ~ _ ‘ 

\±j_  +  irp'"' 


CD  =  (D^  +  ICD^ 


co^  and  cd^  being  real,  and 


1  R2[J  (K  r) ]2 
o v  o  ' 


X2  “  J2. 


R  [j  (K  r)  ] 
o  v  o 


2 


13  R2[Jo(Kor)]2 


'4  R2tJD(V)]2 


R 

/  p  J  (K  r)  dr 
0  1  0 

R  p  p 

I  P  Ji(K0r)  r  dr 


0 


/  |-[r  ^2^  J.  (K  r)  ]  J,  (K  r) 
^  dr L  a  1 v  o  '  1 v  o  ' 

R  2 

I  J{(K0r)  J1(K0r )  r  dr 


To  have  a  non-trivial  solution  for  f  and  &  we 

o  o 

must  set  the  determinant  of  the  coefficients  of  equation 
(22)  equal  to  zero.  That  is 

A ,  B  -  A  B,  =  0  (24) 

~r  ~  “  ~r 

Dividing  equation  (24)  into  its  real  and  imaginary  parts 
we  obtain  two  equations  of  the  form 


Af  cd1  +  cd2  =  Cx 
A  2  ^2  ^2  =  d  2 


(25) 


dr 


. 
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From  these  equations  the  imaginary  frequency  may  be  written 


A 

A, 


1 


C 

C, 


1 


A.  B 


A 


‘1 

2 


B 


At  the  critical  field  the  helical  state  must  neither  grow 
nor  decay.  That  is,  at  the  critical  field,  cd2  =  0.  There¬ 
fore  we  put 

A1  C2  -  A2  C1  =  0  (26) 

recognizing  that  henceforth  Bqz  represents  the  critical 
magnetic  field. 

From  equation  (25) 
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Replacing  k^s  C^,  C2  in  equation  (26),  after 

much  algebra  we  obtain  the  equation 

G^  +  G^  z  +  Gg  z  +  G^  +  [  Gg  z  +  Grj  z  +  Gg  z  ]  ijy 

=  -[  a2  (1  +  y)2  +  4y  ]  t>31  [  G12  (1  +  Cl;L)  z  +  Gl2  z 

+  (  Gg  z  +  G1q  z2  +  Gi;l)  f  ]  (28) 
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and  where  the  coefficients  G  are 
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Setting  if/  equal  to  zero  in  equation  (28),  we  ob¬ 
tain  Johnson's  corresponding  equation. 

For  the  numerical  calculations,  we  shall  consider 
helium  for  which 

P 

—  =  43-7 

q+ 

We  are  now  left  with  an  algebraic  equation  in  the  three 
variables  y,  z  and 

We  wish  to  eliminate  states  which  are  stable  only 
for  a  precise  wavelength  in  favour  of  those  stable  in  a 
small  region  about  the  stable  wavelength.  The  latter  are 
observed  experimentally  but  the  former,  corresponding  to 
a  perfectly  monochromatic  state,  are  not. 

To  impose  this  condition  mathematically,  we  put 
d(lm  a>)  =  n 

or  in  terms  of  the  coefficients, 

6gi  x ^  x^  +  4G2  x^  +  x2  +  2G^x  +  Gg^ 

=  -  [a 2  (1+y)2  +  4yF  b31  [G12(l+Cl;L+  3x2)  +  ‘iG^f  x3 

+  2Ggf  x]  (29) 

1 

where  x  =  z2 . 

Combining  equations  (28)  and  (29)  to  eliminate 
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b^,  which  is  the  only  one  of  the  three  variables  appear¬ 
ing  in  a  simple  enough  form  to  be  eliminated  analytically, 
we  obtain 

2G1G10^  x9  +  3GiGi2  *8  +  (^G!G9  +  2GqGi2)^  x7  +  [5G1G12- 
( lH- C G2G12^  t^GgG]_2  ^"*”G11  ^  GG^G2  "  2G^G-^q 

-6G1G11]^/  x5  +  [3G2Gl2(l+C11)  -  G3G12]  X4  +  [2G7G12* 

(1+C11)  -  2G^G12  -  4G4G10  +  4G2G1X]^  x3  +  [a^G^il+C^) 

-  3G4G12]  x2  +  C2G3Gi;l  -  2G9G4]^/  x  -  G4G12(1+Cll)  =  0  (30) 

Equation  (30)  is  a  relationship  between  z  and  y, 
which  are  essentially  the  wave  number  and  the  critical  mag¬ 
netic  field.  To  evaluate  the  integrals,  p(r)  is  written 
in  terms  of  the  solution  for  n(r)  in  the  steady  state  tak¬ 
ing  only  the  lowest  order  term.  That  is,  we  take 

n(r)  =  J  (y  r) 

Similarly,  EQr  is  replaced  by  its  lowest  order  expression 
from  the  steady  state. 

Dividing  the  range  of  integration  into  ten  equal 
parts  we  obtain  the  following  numerical  values  for  the  int¬ 
egrals  : 

q  =  -  4.307 
i2  =  1.849 
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13  =  -  9.570 

14  =  -  1.625 

From  the  definition  of  ip  we  write 
X 

f  =  (^-)2  S  e  (vz+  -  VZf 

2  ' 


The  orders  of  magnitude  of  the  quantities  involved  are 
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1010  7 

1U  cm/sec 

K 

0 

/v» 

1  /cm 

( Brown 

give 

for 

the  order  of  magnitude 

of  ip 

ip 

IQ"3 

To  accentuate  the  effect  of  the  azimuthal  field, 
equation  (30)  is  plotted  in  figure  2  for  parameter  ip  equal 
to  0.001.  Johnson's  solution,  corresponding  to  parameter 
ip  equal  to  zero,  is  also  plotted  in  the  same  figure. 
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To  obtain  the  graph  of  the  electric  field  as  a 
function  of  the  critical  magnetic  field,  we  must  eliminate 
the  variable  z  between  equations  (28)  and  (29).  This  can 
be  done  only  numerically  by  using  the  points  of  the  graph 
in  x  and  y.  The  relationship  between  the  electric  field 
variable  b^  and  the  critical  magnetic  field  is  plotted 
in  figure  3*  again  for  parameter  f  equal  to  0.001  and  zero. 
This  graph  cannot  be  described  in  an  algebraic  form  from 
equations  (28)  and  (29) . 

The  final  quantity  which  we  wish  to  plot  as  a  func¬ 
tion  of  the  critical  magnetic  field  is  the  frequency  of 
rotation  of  the  helix.  We  must  return  to  equations  (25) 
and  solve  for  cjd^  which  may  be  written 


B 

B 

B 

B 


1 

2 

1 

2 


where  A-^,  A2,  C-^,  are  given  in  (27)  and  where  B1  and  B 2 
are  given  by 


Therefore  may  be  written 
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C-^A-^  +  C^h-2. 

2  ~ 

A1  +  A2 


Substituting  the  values  of  A^,  A^,  Cjy  C2  from  equation 
(27) j  after  much  algebra,  we  obtain  the  following  equation 
relating  the  real  frequency  co^,  the  electric  field  var¬ 
iable  b^2  and  the  wave  number  variable  x  or  z.  This  may 
be  written  as  follows : 
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where  b-^  is  defined  as  follows 


b 


11 


e  a) 

(kT+  +  kTj 


+  4! 


1 

K 


2 

"o 


Using  the  graphical  results  relating  y  and  z  we  calculate 
bn  point  by  point  from  equation  (31).  The  graph  of  b-^  as 


a  function  of  y  is  given  in  figure  4*  again  for  ifs  equal  to 


0.001  and  zero. 
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figure.  2 

Graph  of  the  wavelength  A  of  the  helix  as  a  function 
of  the  critical  longitudinal  magnetic  field  for  para¬ 
meter  ip  equal  to  .001  (solid  curve)  and  zero  (Johnson  !s 
solution,,  dashed  curve) 
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figure  3 


oz 


Graph  of  the  electric  field  variable  b^1=  K  /kT  +kT 

as  a  function  of  the  critical  magnetic  field  for  parameter 

equal  to  .001  (solid  curve-)  and  zero  (Johnson's  solution, 
v. 

dashed  curve) 
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figure  4 

Graph  of  the  frequency  of  the  helix  as  a  function 
of  the  critical  magnetic  field  for  f  equal  to  0.001 
(solid  curve)  and  f  equal  to  zero  (Johnson's  solution,, 
dashed  curve) 
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Chapter  4 :  Conclusions 

4.1  Steady  State 

As  was  mentioned  in  Chapter  2,  it  is  easy  to  ex¬ 
plain  qualitatively  figure  1,  which  shows  the  density  pro¬ 
file  across  the  tube.  It  is  well  known  that  the  azimuthal 
magnetic  field  has  a  pinching  effect  on  the  column  and  its 
inclusion  indeed  results  in  a  higher  concentration  of  ions 
and  electrons  near  the  axis  and  a  correspondingly  lower 
concentration  away  from  the  axis  than  in  Johnson's  calcul¬ 
ation  . 

Verification  of  the  effect  of  the  additional  term 
using  existing  experimental  data  is  impossible.  Experi¬ 
ments  to  date  have  been  performed  in  ranges  of  the  vari¬ 
ables  in  which  the  effect  of  B  is  too  small  to  be  ob- 

ocp 

served.  That  is,  they  have  been  performed  using  low  aver¬ 
age  densities  and  low  currents. 

4.2  Helical  mode 

The  graphs  of  figures  2,  3  and  4  show  the  devia¬ 
tion  from  Johnson's  calculation  for  a  value  of  the  dimen¬ 
sionless  azimuthal  field  parameter  f  equal  to  0.001.  Al¬ 
though  an  experiment  in  this  range  is  possible,  it  would 
have  to  be  one  of  unusually  high  current,  strong  electric 
field  and  high  average  density.  No  experiment  has  been 
performed  to  our  knowledge,  which  was  designed  to  isolate 


the  effect  of  the  azimuthal  field.  The  usual  experiments, 
investigating  the  helical  mode,  are  performed  for  values 
of  y  greater  than  those  where  Johnson's  theory  and  the  pres 
ent  modified  theory  disagree. 

A  representative  example  of  such  experiments  would 
be  one  described  in  a  recent  paper  by  M.  Sato.  This  paper 
gives  a  single  point  on  the  graph  of  wavelength  versus 
critical  field  (figure  2) .  At  the  minimum  of  light  inten¬ 
sity  from  the  column,  which  is  assumed  to  occur  at  the 
critical  field,  a  helix  of  wavelength  A  =  4  R  is  observed 
when  the  critical  magnetic  field  parameter  y  =  10  This 

agrees  with  the  asymptotic  behaviour  of  Johnson's  graph. 

-7 

However,  under  the  conditions  of  the  experiment,  f  =  10 
and  the  deviation  due  to  the  induced  field  in  this  case  is 
important  only  for  y  <  10-^.  Furthermore,  the  observation 
has  no  more  than  order-of-magnitude  precision,  and  since 
the  deviation  between  Johnson's  calculation  and  the  pres¬ 
ent  modification  is  small,  the  only  conclusion  to  be  drawn 
from  this  experimental  point  is  that  it  agrees  with  both 
calculations  within  the  limits  of  its  precision. 

It  is  interesting  to  study  the  behaviour  of  John¬ 
son's  calculation  in  the  limit  that  the  longitudinal  mag¬ 
netic  field  goes  to  zero.  In  Johnson's  calculation,  put¬ 
ting  the  longitudinal  magnetic  field  to  zero  means  there 
is  no  longer  any  magnetic  field  in  any  direction.  The  prob 
lem  reduces  to  a  simple  electric  discharge  in  a  longitudin- 


. 


42. 


al  electric  field  with  a  sink  at  the  boundary.  The  prob¬ 
lem  now  has  cylindrical  symmetry  and  consequently  no  helix 
should  exist.  If  a  helical  state  were  imposed  on  the  sys¬ 
tem.,  it  would  quickly  decay  due  to  the  frictional  force 
retarding  the  azimuthal  motion.  There  is  no  force  which 
could  sustain  the  azimuthal  motion.  The  frequency  and 
the  wavelength  A  may  be  Investigated  in  this  limit  and  it 
is  found  that  ^  tends  to  zero  and  A  tends  to  a  finite  val¬ 
ue.  In  other  words,  Johnson's  calculation  admits  a  sta¬ 
tionary  helix  as  solution  for  Bqz  =  0. 

It  must  be  noted  that  we  have  included  the  in¬ 
duced  field  as  a  perturbation  to  Johnson's  solution.  Con¬ 
sequently,  the  range  of  validity  of  the  modified  results 
is  limited  to  cases  of  small  deviations  from  Johnson's 

L 

results.  The  modification  Is  valid  if  the  deviations  in 
the  wavelength,  electric  field,  and  rotational  frequency 
are  much  smaller  than  the  values  of  the  wavelength,  elect¬ 
ric  field,  and  frequency  respectively. 
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